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I. INTRODUCTION 

In quantum mechanics, the state of a physical system 
is usually described by a density operator, which is a 
positive scmidchnite, Hermitian matrix with unit trace. 
For any pair of Hermitian matrices A and B, define the 
Hilbert-Schmidt inner product (A, B) HS = Tr (AB) . The 
space of Hermitian matrices Herm(H d ) on a Hilbert space 
ri d then forms a Euclidean space. If we consider the set 
of density operators T>CH d ) as a subset of Herm("H d ), 
then we can think of d-dimensional quantum states as 
points in a (d 2 — l)-dimensional convex set C C R d 
that is isomorphic to V(% d ). We expect the geomet- 
ric features of the convex set C to reflect properties of 
density operators. For example, the full geometry of C 
is well-known for d = 2; it is a solid 3-dimensional ball 
called the Bloch ball. The spherical boundary of the ball 
corresponds to pure states, where orthogonal states get 
mapped onto antipodal points, and interior points cor- 
respond to mixtures, each of which can be decomposed 
into any convex combination of pure states whose con- 
vex hull contains that point. We can also compare how 
similar any two states p, a are by measuring their Hilbert- 
Schmidt distance Z?hs(p, c) = y Tr (p — a) . Little, how- 
ever, is known of the same convex geometry for quantum 
states in higher dimensions. Much effort has been made 
in uncovering the rich, intricate structure of C for d = 3 
by examining the various 2- and 3-dimensional sections 
obtained from the generalized Bloch representation for 
qutrits [IH1] but many details of its overall structure re- 
main unknown. 

In this paper, we analyze the geometric features of 
qutrits in terms of the probabilities for a special mea- 
surement called a symmetric informationally complete 
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(SIC) measurement [5T-HU]. A SIC measurement maps 
each density operator into a unique probability vector, 
which represents one way of specifying an isomorphism 
of r D('H d ) onto a subset of E d . This particular mapping 
allows us to characterize quantum states as a proper sub- 
set of the probability simplex, where the restriction is 
imposed mainly by a special version of the Born rule. 

We find that the SIC probabilities provide us with a 
novel way of characterizing properties of qutrits, partic- 
ularly with respect to the flat geometry induced by the 
Euclidean metric defined on the simplex. Here we present 
three main results: 

(i) The extreme points representing pure states for 
qutrits are obtained from a simple formula that 
picks out a submanifold of points lying on a certain 
sphere. 

(ii) The probability vectors obtained for any pair of 
SICs are related by a rotation with a very sim- 
ple form, which we construct explicitly. It may be 
worthwhile to note that the rotation does not neces- 
sarily arise from a unitary transformation between 
the SICs involved (except when the SICs are uni- 
tarily equivalent). 

(iii) The boundary points are described using a polar 
equation that gives their radial distances from the 
uniform distribution, which represents the maxi- 
mally mixed state. This is different from existing 
methods that analyze the geometry of qutrits by 
studying the boundary of various 2-dimensional sec- 
tions. 

The motivation for such a study is twofold. Firstly, the 
geometry of quantum states is interesting in its own right 
and a better understanding of it may have important 
repercussions for various applications of quantum infor- 
mation processing. Secondly, Fuchs and Schack advo- 
cate a framework for reformulating quantum mechan- 
ics directly in terms of probabilities without mentioning 



Hilbert space at all [TT]. A better understanding of the 
structure of SIC probabilities may prove useful in iden- 
tifying the basic axioms needed for reconstructing quan- 
tum theory exclusively in terms of probabilities. 



II. PROPERTIES OF WEYL-HEISENBERG 
QUTRIT SICS 

One way to represent quantum states in terms of prob- 
abilities is to express density operators in terms of d 
linearly independent projectors H = IV'XV'I such that 
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(i) 



When each projection is scaled by ^, we get a measure- 
ment called a symmetric informationally complete posi- 
tive operator- valued measurement (SIC-POVM), a topic 
of considerable interest in the quantum physics commu- 
nity. In this paper, the set {IIi}f =1 is called a SIC for 
short. 

The one-to-one correspondence between the outcome 
probabilities p(i) of a SIC-POVM and density operators 
is given by [T^] 



P = 



E 



(d + l)p(i) 



n, 



(2) 



Thus, the SIC probability vectors p provide an equiv- 
alent description of quantum states, which we call the 
SIC representation. It follows that we can always choose 
the coordinates of C such that p £ C for all p associated 
with a quantum state according to Eq. pi) . The relation 
between the Hilbert-Schmidt inner product for a pair of 
density operators p and a, and the scalar product of their 
respective probability vectors p and q is given by 



Tr (pa) = d(d + l)p-q- 1. 



(3) 



All SICs constructed to date have a certain group co- 
variance property. Let G be a group of d 2 elements and 
let g H> U g be a projective representation of G on H d . 
Let \if>) £ H d . If the set of vectors U g \ip) generates a SIC 
S, we say that S is covariant with respect to G. The seed 
vector \ip) for S is called a fiducial vector. 

In almost all known cases, the unitaries that produce 
SICs belong to the Weyl-Heisenberg group. In d = 3, the 
Weyl-Heisenberg group is generated by the shift X and 
phase Z operators, 




Z = 




(4) 



where ui = e 1 $ . We can act with powers of X and Z on 
a SIC fiducial \i/j) so that the resulting vectors 



\^ mn ) = X m Z n \1>) 



0,1,2, 



(5) 



satisfy Eq. p}. In that case, the projectors associated 
with \rpmn) forms a SIC, which we call a Weyl-Heisenberg 
SIC. 

The operators X and Z generate the group 



W(d) = {uj a X m Z n \ a,m,n= 0, 1, 



1}, (6) 



which we call the Weyl-Heisenberg group. Note that 
W(d) is order d 3 . However, two unitaries which differ 
only by a phase generate the same SIC projector, so the 
SIC itself only contains d 2 elements. As a matter of con- 
vention, we label the Weyl-Heisenberg SIC projectors n^ 
with index i = dm + n + 1, so we have i — 1, 2, . . . , 9 
for qutrit SICs. For example, n 6 is the SIC projector 
corresponding to |*0e) ~ XZ 2 \ijj). 

The normalizer of the Weyl-Heisenberg group is called 
the Clifford group C(d) [8| fT3Hl5] . which is itself a uni- 
tary subgroup in dimension d. If U € C(d) is a Clifford 
unitary operator and W(d) is the Weyl-Heisenberg group 
then 



UW(d)U ] = W(d). 



(7) 



If the set of antiunitary operators that map W(d) to itself 
are included, we get the extended Clifford group. 

Every Weyl-Heisenberg SIC in d — 3 can be obtained 
by acting with an (extended) Clifford (anti) unitary on a 
SIC with fiducial vector 



IV* 




is 



0. 



(8) 



Fiducials corresponding to distinct values of t in the 
range [0, ?] generate distinct orbits of the extended Clif- 
ford group. In Ref. [8 it is shown that there are 3 types 
of orbits of the extended Clifford group in d — 3 for which 
\ipt) is in the orbit: the infinitely many generic ones for 
t £ (0, 5) and 2 exceptional ones for the endpoints t = 
and t = ? . 

6 

In the generic case, each extended Clifford orbit con- 
sists of 8 SICs generated by the fiducial vectors: 



1>: 



(0±) 



m 



(r±) 




(9) 



where r = 1, 2, 3. 

For i = f , there are 4 distinct SICs whose fiducials can 
be chosen as 



v> 



(o) 



v> 



(2) 




v> 



(1) 



V- 



(3) 




(10) 




For t = 0, the fiducial generating the unique SIC can 
be chosen as 



(11) 



It is worth mentioning here that the SICs of Eq. Q are 
inequivalent with respect to Clifford unitaries; however, 
some of the SICs for different values of t are still related 
by a unitary operator that is not a member of the Clifford 
group. Specifically, Zhu [T5] has shown that the SICs 
for t, ? — t and ? + t are, in fact, unitarily equivalent 
to each other, with the unitary transformation relating 
them being 



U = diag(l,u, u ), 



(12) 



which is not a Clifford unitary. Moreover, there are no 
other unitary equivalences. This means that every pair 
of SICs on any two different orbits corresponding to t G 
[0, -jg] are not equivalent. 

Associated with each value of t are two sets of closely 
related geometric quantities. The first set consists of the 
traces of the product of three SIC projectors called triple 
products Tijk, 



T ijk = Tr (11,11,11* 



(13) 



It is show in Ref. [17] that two SICs are unitarily equiva- 
lent if and only if the triple products are the same, up to 
permutation. The other set consists of the structure co- 
efficients Sijk that describe multiplication between SIC 
projectors, 



Urif — 2_^ SijkHk 



(14) 



It is straightforward to show that the structure coeffi- 
cients can be obtained from the triple products in the 
following way: 



^ijk 



1 



K ' tjk d+1 



(15) 



In Sec. Ill we shall see that the real parts Ty* = 
Re [T^k] and Sijk — R- e [Sijk] °f the triple products and 
structure coefficients, respectively, are adequate for de- 
scribing probability vectors corresponding to quantum 
states. 

It is easy to compute Sijk when some of the indices are 
identical: 



&iii — J-) ^ijj — *-*jij — 



1 



Sjji = 0. 



(16) 



It is also straightforward to compute S^k for i =/= j =/= k 
using Eq. (13 1 and Eq. (15). Taking only the real parts, 



the distinct nonzero values are — 4 and 3 other values we 
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TABLE I. Structure coefficient index generators for qutrit 
SICs. The rules for choosing the index triples (ijk) for each 
distinct value of Sijk are described in the main text. 
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(17) 



It can be seen that SICs with parameter values t , ? — t, 



and 



t have the same values of x, y, and 2, up to a 



permutation — confirming the fact mentioned earlier, that 
such SICs arc unitarily equivalent. 

Hughston [IS] has shown that the SIC vectors of the 
single SIC for t = can be obtained from the inflection 
points of a family of cubic elliptic curves on the com- 
plex projective plane known as the Hesse pencil (see also 
Bengtsson |19|.) There are 8 SICs with parameter value 
t — I that are unitarily equivalent to the single SIC for 
t = 0, and we call these 9 SICs the Hesse SICs. The 
particular SIC specified by Eq. (11) shall be called the 
canonical Hesse SIC. 



For the canonical Hesse SIC, Eq. (17) gives 



xo 



ya = z Q = 0. 



(18) 



It is the simplicity of these numbers that leads to an 



elegant characterization of qutrit pure states in Sec. Ill 



There are some simple rules for finding the index triples 
corresponding to the values — \,XfiVu ancl z ti which we 



now describe. 



t 



(r±) 



(r 



XtiVu 
To each qutrit SIC-family 

0, 1, 2, 3) in Eq. ([£]) we assign an index generator G r ± 
that helps us choose the index triples (ijk) for each dis- 
tinct value of Sijk' They are listed in Table [El To demon- 
strate what the rules are, it is helpful to consider the SIC 

ip\ ) as a specific example. 

For S^k — — i) take the index triples on the same row. 
Looking at G2+ in Table III we see that the relevant set 



of {ijk) for 



t. 



(2+) 



{(168), (249), (357))} 



and all permutations of indices for each {ijk). 

For Sijk = x t , take the index triples belonging to the 
same column, or those on entirely different rows and 
columns. Thus, the relevant set of (ijk) from G2+ is 

{(123), (645), (897), (147), (693), (825), (195), (627), (843)} 

and all permutations of indices for each (ijk). 

For Sijk = Vu take the index triples such that first two 
indices belong to the same column, and the last one is 
in a different row and belongs to the succeeding column 
when counting in a cyclic manner. By succeeding we 
mean that "column 2 is after column 1 " , "column 3 is 
after column 2" , and "column 1 is after column 3." Thus, 
the relevant set of (ijk) from G2+ is 

{(125), (647), (893), (236), (458), (971), (314), (569), (782)} 

and all permutations of indices for each (ijk). 

For Sijk — zt, we have a similar rule as in y t but take 
the last index from the preceding row. Thus, the relevant 
set (ijk) from G2+ is 

{(127), (643), (895), (238), (451), (976), (319), (562), (784)} 

and all permutations of indices for each (ijk). 

Any other index triple (ijk) not specified above has 

Sijk = U. 

We do not need them in this paper but it is possible to 
construct similar, though somewhat more complicated, 
rules for getting the imaginary parts of Sijk- 



III. PURE STATES IN THE SIC 
REPRESENTATION 

Since quantum state space is a compact convex body 
C in M. d , the Krcin-Milman theorem 20 1 states that it is 
equal to the convex hull of its extreme points, which are 
the pure states. It is therefore natural to ask what are 
the conditions on SIC probability vectors p such that they 
correspond to pure states. In terms of density operators, 
a pure state is represented by a rank-1 projector, p 2 = p. 
A remarkable theorem [5TJ [53] states that for a hermitian 
operator p — p\ an equivalent condition for defining a 



pure state is given by 



Tr(p 2 )=Tr(p 3 )=l. 



(19) 



Using the SIC representation of p given by Eq. (pi), 
Eq. Jl9| becomes [231 



i 

E T ljk p(i)p(j)p(k) 



i,j,k 



d(d+l)' 

d + 7 
(d+1) 3 ' 



(20) 
(21) 




FIG. 1 . The 12 lines of a finite affinc plane over the Galois field 
GF(3) representing the index triples (ijk) £ Q in Eq.(27l. 
The indices are depicted as 9 points and each line contains 3 
points as marked. 



Since the right-hand side of Eq. (21) is real, and since 
the imaginary parts of the triple products are completely 
antisymmetric, we have 

E T ljkP (i)p(j)p(k) = ]T fijk P (i)p(j)p(k), (22) 



i,j,k 



i,j,k 



so we may consider just the real parts Tijk- We obtain an 
equivalent expression for Eq. (21) in terms of the struc- 
ture coefficients: 



E Sijkp(i)p(j)p(k) = 



i,j,k 



d(d+l) 2 



(23) 



Specializing to the case d = 3, the pure states are 
described by probability vectors p that satisfy 



1 



i 

E S ljk p(i)p(j)p(k) = — . 



i,3,k 



(24) 
(25) 



Using the results in Sec. |TTJ we find that the pure states 
for the canonical Hesse SIC are given by 

1 
6' 



Epw 2 



(26) 



E^) 3 = 9 S p( i )p(tip( k ) 



(27) 



(ijk)eQ 



where Q is the set of index triples (ijk) corresponding 
to the lines drawn in Fig. [IJ where permutations of the 
indices i,j,k are counted separately. Interestingly, Fig. [T] 
coincides with a combinatorial object called a finite affine 
plane of order 3 (which can also be identified with the 
unique 2-(9, 3, l)-design) [24]. It contains 9 points and 
12 lines, and the index triples in Q correspond to any 3 
points on the same line. 

The pure states for any other qutrit SIC are located on 
the same sphere given by Eq. J24"l ) but with a different 
set of values for Sijk m E q. (|23f 7 If we substitute the 
values in Eq. ([l6]) into Eq. (|23|, we obtain 



oE*(* 



E Si]kP(i)p(j)p(k) 



0. 



(28) 



IV. ORTHOGONAL TRANSFORMATIONS 

BETWEEN PROBABILITY VECTORS OF 

DISTINCT SICS 



In this section, we will show that, in any dimension, the 
probability vectors corresponding to two different SICs 
are related by an orthogonal transformation. We will 
then go on to show that in dimension 3, the orthogonal 
transformation takes a remarkably simple form. 

In arbitrary dimension d, consider a pair of distinct 
SICs with elements U^ and ILy . In the vector space of op- 
erators, 11^ and EL, correspond to the vertices of 2 identi- 
cal regular simplices, which means there must be related 
by an orthogonal transformation. Formally, because a 
SIC forms a Hermitian basis in the space of operators, 
we can write 



ii^E^il, 



(29) 



Since every SIC element has unit trace, taking the trace 
on both sides of Eq. ( 29 ) gives us 



E^ = L 



(30) 



Multiplying the left-hand side of Eq. ( 29 1 by IK, and the 
right-hand side by J2i Rlj^-it we nave 



fdS k i + 1 



Tr(n;n;)=^ Wj! 

k,l ^ 

= '* ^i = / ,RikRjk- 



(31) 



which confirms that 7?j, is indeed an orthogonal matrix. 
Using Eq. (31), it is now straightforward to show that 



p'(i) = J2RnP(i)- 



(32) 



Let Hi be the canonical Hesse SIC and let III (t) be 



(r±), 



the SIC generated by the fiducial vector tf>. 



= 3 (4n. t - 1) 



dr±) 



be the dual basis to IE. (so Tr(QiUj) — 5 Ki 
orthogonal matrix which takes IE onto IE; (t) is 



R^(t) = Tr IL^(t)Qj 



Let 

(33) 
Then the 



(34) 



It turns out that the matrices R( r± ' (t) have a very simple 



form. Define the permutations 
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Let P^ r± > be the permutation matrix corresponding to 
p( r± > , with matrix elements 



>(r±) _ 



3j,p(r±) (i ). 



Also define 



o(t) = -(l + 2cos2i), 



(36) 



(37) 



a(t) a(t- f) a(t + f) 
A(t) = | a (t + |) o(*) a (t - f ) | , (38) 



l0 ft-f) o(t+f) 



o(t) 



^(t) 
i?(t) = | A(t) 
A(t) 



(39) 



It is then straightforward, though somewhat tedious, to 
verify that 



R ir± \t) 



p(r±) 



R{t)P (r±) 



(40) 



Since R(t) = I ® A(i), it follows that 

Det[i? (r±) (i)] = Det[i?(t)] = {Bct[A(t)]} 3 . (41) 

Because A(t) is a circulant matrix, its eigenvalues are 
given by 

A m = a(t) + cj m a (t - |) + w- m a (t + |) 

= e 2ttm (42) 

for m = -1, 0, 1. This implies that Det[A(*)] = 1. Thus, 
Det[i?( r± )(i)] = 1 and i?( r± )(i) is, in fact, a rotation 
matrix. 

It is easily seen that 

R(h)R(t 2 ) = R{t 1 + t 2 ), 72(0) = 7. (43) 

So the matrices R(t) form a 1-parameter subgroup of the 
orthogonal group. 



V. THE BOUNDARY OF QUTRIT STATE 
SPACE 

A concrete way to understand the geometry of qutrit 
state space is to figure out what the convex body looks 
like. In this regard, we want to consider not just the 
pure states but all boundary points of the set. Some 
valuable insight into the shape of the boundary is gained 
by looking at the distance of the boundary states from the 
center of the space, the maximally mixed state p = \l 
as a function of direction. Specifically, we can write the 
SIC probabilities in the form 



P(i) = -ip +rn(i) 
where ft is a direction vector with 



5>(*) = o, 5>(i) 2 = i. 



(44) 



(45) 



Let r(ft) be the value of the r corresponding to the quan- 
tum state on the boundary, which is given by Eq. ( 44 1 



Here we calculate this function for the canonical Hesse 
SIC. Put differently, we are looking for the polar equation 
describing its boundary states. 

The boundary is determined using the following 
lemma: 

Lemma 1. Let p be an arbitrary Hermitian operator on 
a 3 -dimensional Hilbert space. Then 

(i) p is a density operator if and only if 

Tr (p) = 1, Tr (p 2 ) < 1, 3Tr (p 2 ) - 2Tr (p 3 ) < 1; 

(ii) p is a density operator for a boundary state if and 
only if 

Tr (p) = 1, Tr (p 2 ) < 1, 3Tr (p 2 ) - 2Tr (p 3 ) = 1; 

(Hi) p is a density operator for a pure quantum state if 
and only if 

Tr (p) = 1, Tr (p 2 ) = 1, 3Tr (p 2 ) - 2Tr (p 3 ) = 1. 

Proof. We begin by proving necessity. Suppose p is a 
density matrix. It immediately follows that Tr(p) = 1 
and Tr(p 2 ) < 1. To prove the remaining inequality, let 
a, j3, 1 — a — p be the eigenvalues of p. We find 

3Tr(p 2 ) - 2Tr(p 3 ) - 1 = -6af3(l -a-f3)<0. (46) 

For p to be a boundary state at least one of its eigenvalues 
must vanish, in which case 

3Tr(p 2 ) - 2Tr(p 3 ) -1 = 0. (47) 

In addition, if p is a pure state then Tr(p 2 ) = 1 



We now turn to the proof of sufficiency. Let p = p^ be 
such that 

Tr(p) = l, Tr(p 2 )<l, 3Tr(p 2 ) - 2Tr(p 3 ) < 1. (48) 

The first equality means that we can take the eigenvalues 
oi p to be a, (3,1 — a — j3. From Eq. (46 ) we get 



0^(1- a- £) >0. 



(49) 



Thus, either (i) all eigenvalues are non-negative or (ii) 
exactly two of them are negative. We can show that (ii) 
is impossible. Assume the contrary to hold. Without 
loss of generality a, /3 < 0, implying that 1 — a — fj > 1, 
which in turn implies Tr(p 2 ) > 1, contrary to hypoth- 
esis. We conclude that p is positive semi-definite, and 
consequently a density matrix. 
Next assume that 

Tr(p) = 1, Tr(p 2 ) < 1, 3Tr(p 2 ) - 2Tr(p 3 ) = 1. (50) 

Then, 

a(3{\ - a - (3) = 0, (51) 

implying that at least one of the eigenvalues must be 
zero. So p is on the boundary of state space. 
Finally assume 

Tr(p) = l, Tr(p 2 ) = l, 3Tr(p 2 ) - 2Tr(p 3 ) = 1. (52) 

Using the argument above the eigenvalues are 0, a, 1 — a. 
Since Tr(p 2 ) = 1 it must be that a = or 1 and therefore 
p is a rank-1 projection operator. □ 

We can use the lemma for the quantum states associ- 
ated with the canonical Hesse SIC. To this end, recall its 



structure coefficients in Eq. (18). We use these to calcu- 
late Tr (p 2 ) and Tr (p 3 ) for p given by Eq. J2~l). We find 
that 



Tr(p 2 )=12]Tp(,) 2 -l, 



(53) 



Tr(p 3 )=l + 24j>(z) 3 -12 £ p(i)p(j)p{k) , 

(ijk)eQ 

where Q is again the set of lines on the affine plane in 
Fig.[l| Substituting Eq. (44 1 into the probabilities above, 
we obtain 



i 
i 



(ijk)£Q 



r 

y 

r 2 
3 



■E-» 3 



(54) 



r y n(i)n(j)n(k) . 

(ijk)£Q 



Consequently, we can restate the conditions in Lemma [T] 
as follows: 




FIG. 2. A plot of F as a function of one of the eigenvalues a 
for boundary states. 



(i) p is a density operator if and only if 



r z < 



18' 



Ar 3 F(n) -r 2 -\ > 0; 

54 ~~ 



(55) 



(ii) p is a density operator on the boundary of the state 
space if and only if 



r 2 < 



i 4r 3 F(n)-r 2 + ~=0; (56) 



(iii) p is a density operator for a pure state if and only 
if 



18' v ' 



1 



where 



F{n) = J]n(i) 3 - - 53 n(i)n{j)n(k). 



(57) 



(58) 



(ijk)ec 



Thus, the value of r(n) giving the distance of a boundary 
state from the completely mixed state along the direction 
of n is the smallest positive root of 



Ar 3 F(n) - r 2 H = 0. 

54 



(59) 



To determine the bounds on F = -F(n), we can write it 
in terms of the eigenvalues of p by performing some alge- 
bra on Eq. (53) and Eq. (54). Denoting the eigenvalues 



of p by a, (3, and 1 — a — p, we obtain 



n/3| 



/1 + /2 



(/1 " I)"* 

where < a, /3, 1 — a — /3 < 1 and 

h=a 2 + p 2 + (l-a-P) 2 , 
f 2 =a 3 + f3 3 + (l-a-f3) 3 . 



(60) 



From Eq. ( 60 ) , it can easily be shown that 
-±<F<±- 

V2~ ~ y/2 



(61) 



(62) 
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FIG. 3. The radial distance r of boundary states from the 
maximally mixed state, as a function of F = F(n). 



where the upper (respectively, lower) bound is achieved 
when two of the eigenvalues are identical and < | (re- 
spectively, > I ) . If all the eigenvalues are equal to | this 
corresponds to the maximally mixed state, for which F is 
undefined. For boundary states, at least one of the eigen- 
values must be zero. So the only case we need to consider 
for F is when (3 = 1 — a. Fig. [2] shows F as a function 
of a, provided that one of the eigenvalues vanishes. 



In terms of F, the desired root in Eq. ( 59 ) is given by 



-i-(l+ w - s . g +-) ifF^O 
r = < U f V 9 J (63) 

37! 

where s — sgn(i^) and g is the cube root with the smallest 
positive argument in 



g 3 = 1 - AF 2 + 2F^AF 2 - 2. 



(64) 



When F — ±-75, g = e 1 ^ and we get the bounds for r: 

(65) 
A plot of r as a function of F is shown in Fig. [3] 



1 1 

6\/2 ~ r ~ 3a/2" 



VI. SUMMARY AND OUTLOOK 

With the recent revival of interest in addressing foun- 
dational issues in quantum theory, we pose a simple yet 

intriguing question: what is the shape of the set of quan- 

,2 
turn states C C H. ? Our preliminary attempt to address 

this question revolves around a description of quantum 
states in terms of the outcome probabilities of a SIC- 
POVM. This particular representation allows us to ex- 
ploit the intrinsic symmetry of a SIC in mapping den- 
sity operators to probability vectors, which we believe 
not only serves a natural "coordinate system" for study- 
ing the underlying geometry of quantum states, but also 



provides us with an interpretation of quantum states in 
terms of Bayesian probabilities [25] . 

In this work, we focused our attention on d — 3, which 
is the simplest, nontrivial case to examine. We consid- 
ered the infinitely many Weyl-Heisenberg qutrit SICs, 
which are classified into SIC-families corresponding to 
orbits of the Clifford group. Each SIC can be uniquely 
identified with a set of complex numbers called triple 
products Tijk, the trace of the product of three SIC el- 
ements, whose polar angles are related to discrete geo- 
metric phases J2H1 H7] and to Bargmann invariants [25] 
in complex projective space, and whose imaginary parts 
give SICs the structure of a Lie algebra [TTj. We also 
have structure coefficients Sijk, which are the expansion 
coefficients when multiplying SIC projectors, and whose 
real parts Sijk are especially convenient for describing 
geometric properties of qutrits. 

Using S^k for the canonical Hesse SIC given by \ipo) 
in Eq. (11), we discovered the most economical descrip- 



tion for SIC probability vectors associated with qutrit 
pure states, which are given by Eq. (26 1 and Eq. (27 1. 



Studying the probabilities for the canonical Hesse SIC 
is sufficient because we demonstrated that the proba- 
bilities for other qutrit SICs are related to it by a 9- 
dimensional rotation that can be expressed in terms of a 
single function — a(t) in Eq. d37b. 



The remarkable simplicity of Eq. ( 27 ) suggests that the 



geometric structure of qutrits is largely determined by 
the symmetries associated with a finite affine plane. For 
example, observe that if we consider the indices as points 
in Fig. [II the permutations given in Eq. ( 35 ) is such that 
they preserve the affine lines. Also, using the notion of 



maximal consistent sets in Ref. [39] , a set of probability 
vectors associated with the Hesse configuration of vectors 
on a Hilbert space can be maximized into a convex body 
that has the same largest inscribed sphere and smallest 
containing sphere as qutrit state space, and also shares 
some of its 2-dimensional sections. Therefore, to gain a 
proper understanding of the convex geometry of qutrits, 
it is crucial to understand the full significance of Eq. ( 27 1 . 



We also described a polar equation for the qutrit 
boundary. We found that the function F(n) that yields 
the radial distance of a boundary state in direction n 
from the uniform distribution is the same function that 
picks out the pure states for the canonical Hesse SIC, i.e., 
F(p) = for any pon the sphere containing pure states, 
which again highlights the important role played by the 
finite affine plane of Fig. [l] 

Lastly, it is our hope that the results presented here 
also serve as evidence for the utility of the SIC represen- 
tation in matters regarding quantum foundations. 
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